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Potential Models for Quarkonia 



What can we learn from potential models at finite temperature? 
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Abstract. In this paper I discuss what we can learn about quarkonium dissociation from lattice-potential 
based models. Special emphasis is given to results obtained in agreement by different models, and to the 
relevance of lattice QCD for potential models. Future directions are also discussed. 
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Why are we interested? 

One of the aims of relativistic heavy ion collisions is to 
produce quark-gluon plasma (QGP), a state of matter 
in which the constituents of our hadronic world are de- 
confined. Deconfinement is expected to happen at large 
energy densities, which can be obtained by heating mat- 
ter to extreme high temperatures. High energy density 
matter has been already produced at SPS CERN and at 
RHIC BNL, and will be produced at the LHC, which just 
started its operation at CERN. To have control over what 
temperatures are achieved and whether deconfined matter 
has been produced we need a thermometer. The sequen- 
tial melting of quarkonium has been long considered to be 
exactly that: the QGP thermometer [I]. In a deconfined 
matter the force between the constituents of a quarkonium 
state, a heavy quark and its antiquark, is weakened by the 
color screening produced by the light quarks and gluons. 
For twenty years it has been believed that this screen- 
ing leads to the dissociation (melting) of quarkonium [2]. 
The different quarkonium states are expected to melt se- 
quentially, at different temperatures. A suppressed yield 
of quarkonium can be visible in the dilepton spectrum, 
which is measured in experiments. 

J/tp suppression has been indeed measured by the dif- 
ferent experiments [3J. Understanding the data, however, 
turned out to be more complicated. The reason is that 
the suppression pattern seen is not only due to the hot 
medium effects of screening, but more like due to the in- 
terplay of this with effects of cold nuclear matter [4], as 
well as those of recombination 5 . In order to disentan- 
gle these different effects we must know the properties of 
quarkonium in-medium and determine their dissociation 
temperatures. 
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In principle, everything about a given quarkonium chan- 
nel is embedded in its spectral function: The position of a 
peak in the spectral function corresponds to the mass of a 
bound state, while its width determines its lifetime. Melt- 
ing of a state corresponds to the disappearance of a peak. 
A spectral function also contains information about the 
continuum and its threshold. So following how the spectral 
function changes with temperature can give us a theoreti- 
cal insight to the temperature-dependence of quarkonium 
properties. There are two main lines of theoretical stud- 
ies to determine quarkonium spectral functions at finite 
temperature: potential models and lattice QCD. Poten- 
tial models have been widely used to study quarkonium, 
but their applicability at finite temperature is still under 
scrutiny. Lattice QCD provides the most straightforward 
way to determine spectral functions, but the results suffer 
from discretization effects and statistical errors, and thus 
are still inconclusive. 

In the rest of the paper I discuss our current under- 
standing of what we can learn using the potential models 
at finite temperature currently on the market, highlight- 
ing the agreement produced by the different groups, as 
well as the usage of lattice data (as input or as constrain) 
in context of these models. 



What Goes Into Potential Models 

Potential models are based on the assumption that the 
interaction between the heavy quark and its antiquark in- 
side the quarkonium can be described by a potential. Due 
to the largeness of the heavy quark mass, m c ,& ^> Aqcd 
and the smallness of the heavy quark velocity, one 
can treat the quark-antiquark system nonrelativistically 
and solve the Scrodinger equation to obtain the bound 
state properties. 
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The zero temperature potential model with the Cornell 
potential has experienced great success: It describes well 
the experimentally observed quarkonium spectroscopy [6] ; 
It is verified on the lattice [7]; and it can be derived di- 
rectly from QCD [5]. The latest is possible due to the hi- 
erarchy of well separated energy-scales m 3> mv ^> mv , 
which allows to systematically integrate out the different 
scales and obtain the non-relativistic potential QCD (pN- 
RQCD) (in this the Cornell potential shows up as the 
zeroth order matching coefficient) [8]. 

Inspired by its success at zero temperature the poten- 
tial model has been applied at finite temperature, with the 
main assumption that medium effects can be accounted for 
as a temperature-dependent potential. A few years ago it 
has become clear that instead of just looking at the in- 
dividual bound states (procedure good at T = where 
quarkonium is well defined), we should rather obtain a uni- 
fied treatment of bound states, threshold and continuum 
by determining the spectral function [9j. Spectral func- 
tions can be obtained in a number of ways: by determining 
all the discrete states from solving the Schrodinger equa- 
tion [10] : or by using a T-matrix approach [IT] : or, and this 
has been our choice, by solving the Schrodinger equation 
for the nonrelativistic Green's function, and then make 
use of the optical theorem [12]. In all these approaches 
one has to correctly account for the relativistic part of the 
continuum |12j . 

At first, potential models at finite temperature used 
a phenomenological screened Cornell potential [13]. This 
is all right for merely qualitative purposes. In order to 
get quantitative about how the bound state properties 
(e.g. binding energy) change with temperature, knowledge 
of the exact potential at finite T is required. In princi- 
ple, one should derive this directly from QCD, just as 
this has been done at zero temperature. Such derivation, 
however, is complicated by the existence of temperature- 
driven scales, T, gT, g 2 T , and has only been addressed re- 
cently P3] . So in lack of knowledge resuming to different 
phenomenological potentials came naturally. In particular, 
it became popular using so-called lattice-based potentials, 
i.e. potentials constrained by lattice data on the free en- 
ergy of a static quark-antiquark pair. More precisely, by 
the change in the free energy of a medium at a given tem- 
perature when a static quark-antiquark pair is immersed 
into it. The lattice results show that above deconfinement 
the range of interaction between the quark and antiquark 
is strongly reduced, and this effect can be well described 
by exponential screening [15] . 

One of the most debated questions of recent years is 
which lattice-based potential is to be used in the Schrodinger 
equation. First, it was the free energy F\ [16] . Then, it has 
been understood that this serves merely as a lower limit, 
since it contains an entropy contribution [15] . Removing 
the entropy the internal energy U\ is obtained 15 . U\ 
has been then used as potential |llj . The interpretation 
of U\ as potential is also questionable^, thus serving only 
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Ui includes medium polarization effects; has a huge in- 
crease near T c ; and has an unmotivated increased strength at 
short distances compared to the T — potential [15] . 
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Fig. 1. Set of potentials allowed by lattice data on free ener- 
gies. (See [9] for details) 



as a sort of upper limit. Other lattice-based potentials 
on the market include the one proposed by Wong as a 
combination of F\ and U\ [T7lfT0]. and a set of potentials 
constructed using the general features of the lattice free 
energy: at short distances no deviation form the vacuum 
potential, at large distances exponential screening 12J. A 
set of lattice-based potentials is shown in the left panel of 
figure [TJ 



What Comes Out of Potential Models 

It is clear that different potentials can lead to different 
spectral functions and thus different properties for the 
quarkonium states. More precisely, the deeper the poten- 
tial is the stronger the binding of quarkonium can be, and 
the higher the possible dissociation temperatures might 
then be. 

Granted the differences, there are, however, a number 
of essential features common for all the spectral functions 
from potential models [T0 l fTT [ fT2j : 

1) There is a large threshold (rescattering) enhancement 
beyond what corresponds to free quark propagtion. This 
enhancement is present even at high temperatures, and 
it is the indication that there is correlation persisting be- 
tween the quark and antiquark. Threshold enhancement 
has been identified in all of the channels (charmonium and 
bottomonium S- and P-states). 

2) There is a strong decrease with increasing temperature 
of the binding energies (the distance between peak posi- 
tion and continuum threshold in the spectral function) de- 
termined from potential models. This is illustrated in Fig- 
ure [2J for the J/ip (upper panel) and the T (lower panel). 
The black line on the figure corresponds to the tempera- 
ture line, where the binding energy becomes of the order 
of temperature. 

What we learn from this is that one can obtain a spec- 
tral function that exhibits a resonance-like peak, but the 
corresponding binding energy can be small. It becomes 
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Fig. 2. The binding energy of the J/ip (upper panel) and the 
T (lower panel) determined in different potential models. The 
black line corresponds to the temperature line. 



meaningless, for instance, to talk about a 3 GeV mass 
J/-0 with 20 MeV binding energy! Therefore, unlike for 
quarkonium at zero temperature, at high temperature it is 
incorrect to call a state dissociated only when it's binding 
energy becomes zero. The condition for a state to be con- 
sidered dissociated must be thus weaker than Eu n = . 
We discuss this in more detail later in this paper. At this 
point let me emphasize, that with increasing temperature 
a state can become quickly broadened (weakly bound) and 
thermal fluctuations can dissociate it. Thermal broaden- 
ing of quarkonium, although not taken into account in 
potential model spectral functions, has been addressed in 
a number of independent calculations: NLO perturbative 
QCD [18], QCD sum rule [19], and resummed perturbative 
QCD [20]. All of these calculations show that, for example, 
the J/ip is significantly broadened at temperatures right 
above that of deconfinement. 

It is important to keep in mind when looking at spec- 
tral functions from potential models using lattice-based 
potentials: none of these calculations include the true width 
of a state. Therefore just seeing the peak structure in these 
spectral functions is incomplete on its own and can be mis- 
leading. Note however, that in the perturbative approach 



the potential has also an imaginary part [T4ll20] l2T] . And 
the inclusion of the imaginary part clearly broadens the 
peak [20]. 

Besides illustrating the decrease of binding energies 
with increasing temperature, Figure O also shows that 
there are large uncertainties from the modeling of the po- 
tential. In the following I will discuss whether we can dis- 
tinguish which potential is the "right potential" , if any, by 
making further use of lattice data. 



Potential Models & Lattice QCD 

We have discussed that in lack of knowing the finite tem- 
perature quark-antiquark potential in QCD (pNRQCD) 
the free energy of a static quark-antiquark pair calcu- 
lated on the lattice inspired a series of lattice-based po- 
tentials used as input for potential models. There are, 
however, two independent lattice QCD calculations rel- 
evant for quarkonium studies: One is the calculation of 
the T-dependence of the free energy of a static quark- 
antiquark pair [15] , and the other is the calculation of the 
current-current correlation function of mesonic currents 
in Euclidean-time. From the current-current correlation 
functions for different quarkonium channels the quarko- 
nium spectral functions are extracted [22] . The extraction 
of spectral function using the Maximum Entropy Method 
is still difficult due to discretization effects and statistical 
errors. So the uncertainties in the spectral function are 
significant and details of this cannot be resolved! Thus it 
is difficult to make any conclusive statement based on the 
shape of the spectral functions 23J . 

Contrary to the extracted spectral functions, the lat- 
tice correlators are determined reliably. Therefore, it has 
been suggested to compare correlators from potential mod- 
els to the ones from lattice QCD [5]. Correlators in poten- 
tial models can be calculated using the integral represen- 
tation of these in terms of the spectral function, 

G(t,T) = J du;a{uj,T)cosh(uj(T-l/2T))/smh{uj/2T) . 

It is customary to present the temperature-dependence of 
the correlators as the ratio G/G recon , where G recorl (r, T = 
0) is the correlator reconstructed from spectral function 
at zero (or some low) temperature. This ratio evaluated 
in potential models should then be compared to this ratio 
directly calculated on the lattice. 

The chart in Figure [3] shows a diagram of the inter- 
connectedness of lattice QCD and potential models. Note 
though, that these approaches describe quarkonium spec- 
tral functions in equilibrated QGP. In order to bridge 
these theoretical approaches to experimental data one needs 
to resume to some dynamical model (see for instance [24]). 

So the idea is to constrain different potential models 
using data from lattice [T2]. The approach is the following: 
Assume a lattice-based potential and determine the corre- 
sponding spectral function in a given quarkonium channel. 
From the spectral function determine the ratio of correla- 
tors G/G rec and compare it to lattice data. The "correct" 
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Fig. 3. Structural chart of lattice QCD and potential model calculations. 



potential and spectral function would be the one in best 
agreement with the data. The surprising result is shown in 
Figure [4] This figure displays the correlator ratios in the 
pseudoscalar charmonium channel at 1.2T C obtained us- 
ing the set of potentials within the allowed ranges shown 
in Figure [T] (except for U\), and compared to lattice data 
from [23]. The surprise is that the complete set of po- 
tentials all provide agreement of 1 — 2% accuracy with 
lattice correlators, yielding indistinguishable results. Us- 
ing the internal energy U\ as potential yields worst re- 
sults [TT]. Interestingly, spectral functions that exhibit a 
ground state peak up to higher temperatures, obtained 
with the more binding potentials, and spectral functions 
with no resonance- like structure seen already near T c , ob- 
tained with less binding potentials, (see spectral functions 
in [T^]) all yield flat correlator ratios. This is understood as 
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Fig. 4. The correlator ratio obtained with the set of potentials 
from Fig. [1] compared to lattice data (right) at 1.2T C . (See [9] 
for details) 



the threshold enhancement compensates for the melting of 
states keeping the correlators fiat. Thus possible dramatic 
changes in the spectral function are not reflected in the 
correlator. We therefore cannot identify the "correct" po- 
tential, cannot determine the exact spectral function and 
quarkonium properties from such comparisons. 

The agreement to the correlator ratio is also good for 
bottomonium S-wave, as well as for all the P-states |12j . 
It is now understood that the dominant source of the T- 
dependence of the correlator ratios comes from the com- 
monly overlooked zero- modes [25) . These are low- frequency 
contributions to the spectral functions at finite T describ- 
ing the scattering states of single heavy quarks, in ad- 
dition to the usual bound and unbound quark- antiquark 
pairs in a given channel. The zero-mode contribution is 
understood in terms of quasi-free quarks with some ef- 
fective mass [26], indicating the presence of free heavy 
quarks in the deconfined phase. Furthermore, the zero- 
mode gives a constant contribution to correlator. One can 
eliminate it by looking at ratio of the derivatives of corre- 
lators [25j . The high energy part which carries info about 
bound states is flat, i.e. shows almost no T-dependence 
until 3T C in all of the channels! One can conclude that 
the flatness is not related to survival, since it would also 
imply that the Xc survives until 3T C . The understanding 
becomes then simple: the dramatic changes seen in the 
spectral function are not reflected in the correlator. In 
case of the T this means that the spectral function does 
not change to higher temperatures, the state survives, and 
the correlator ratio stays flat. For all the other channels 
(J/tpi Xc-, Xb-, etc) the spectral enhancement near threshold 
compensates for the melting of states keeping the corre- 
lators flat. Let me make this statement more transparent 
through the discussion in the next section. 

So What Have We Learned? 

As discussed earlier, it has been customary to consider a 
state dissociated when its binding energy becomes zero. 
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Fig. 5. The QGP thermometer. 



In principle, a state is dissociated when no peak struc- 
ture is seen, but the widths shown in spectral functions 
from current potential model calculations are not physi- 
cal. Broadening of states as the temperature increases is 
not included in any of these models. At which T the peak 
structure disappears then? In [57] we argue that no need 
to reach Eu n — to dissociate, but when E^m < T a state 
is weakly bound and thermal fluctuations can destroy it. 
Let us quantify this statement. 

Due to the uncertainty in the potential we cannot de- 
termine the binding energy exactly, but we can never- 
theless set an upper limit for it [27] : We can determine 
Ebin with the most confining potential that is still within 
the allowed ranges by lattice data on free energies. For 
the most confining potential the distance where deviation 
from T = potential starts is pushed to large distances 
so it coincides with the distance where screening sets in 
[12] . From Eun we can then estimate, following [28], the 
quarkonium dissociation rate due to thermal activation, 
obtaining this way the thermal width of a state r(T). 
At temperatures where the width, that is the inverse of 
the decay time, is greater than the binding energy, that is 
the inverse of the binding time, the state will likely to be 
dissociated. In other words, a state would melt before it 
binds. For example, already close to T c the J/V> would melt 
before it would have time to bind. To quantify the dissoci- 
ation condition we have set a more conservative condition 
for dissociation: 2E bm (T) < T(T). The result for differ- 
ent charmonium and bottomonium states is shown in the 
thermometer of figure Note, that all these numbers are 
to be though of as upper limits. 

In summary, potential models utilizing a set of poten- 
tials between the lower and upper limit constrained by 
lattice free energy lattice data yield agreement with lat- 
tice data on correlators in all quarkonium channels. Due 
to this indistinguishability of potentials by the data the 



precise quarkonium properties cannot be determined this 
way, but the upper limit can be estimated. The decrease 
in binding energies with increasing temperature, observed 
in all the potential models on the market, can yield sig- 
nificant broadening, not accounted for in the currently 
shown spectral functions from these models. The upper 
limit estimated using the confining potential predicts that 
all bound states melt by 1.3T C , except the Upsilon, which 
survives until 2T C . The large threshold enhancement above 
free propagation seen in the spectral functions even at high 
temperatures, again observed in all the potential models 
on the market, compensates for melting of states (yielding 
fiat correlators), and indicates that correlation between 
quark and antiquark persists. Lattice results are thus con- 
sistent with quarkonium melting. 



And What's Next? 

Implications of the QGP thermometer of figure[5]for heavy 
ion collisions should be considered by phenomenological 
studies. This can have consequences for the understanding 
of the i?AAUieasurements, since now the J?/> should melt 
at SPS and RHIC energies as well. The thermometer also 
suggests that the T will be suppressed at the LHC, and 
that centrality dependence of this can reveal whether this 
happens already at RHIC. So measurements of the T can 
be an interesting probe of matter at RHIC as well as at 
the LHC. 

The exact determination of quarkonium properties the 
future is in the effective field theories from QCD at finite 
T. First works on this already appeared [T4] and both real 
and imaginary parts of the potential have been derived 
in certain limits. In these works there is indication that 
most likely charmonium states dissolve in QGP due ther- 
mal effects, such as activation to octet states, screening, 
Landau-damping. 

The correlations of heavy-quark pairs that is embedded 
in the threshold enhancement should be taken seriously 
and its consequences, such as possible non-statistical re- 
combination taken into account in dynamic models that 
attempt the interpretation of experimental data [24] , 

All of the above discussion is for an isotropic medium. 
Recently, the effect of anisotropic plasma has been con- 
sidered [29] , Accordingly, quarkonium might be stronger 
bound in an anisotropic medium, especially if it is aligned 
along the anisotropy of the medium (beam direction). 
Qualitative consequences of these are considered in an up- 
coming publication [3D]. Also, all of the above discussion 
refers to quarkonium at rest. Finite momentum calcula- 
tions arc under investigation. It is expected that a moving 
quarkonium dissociates faster. 
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